Anomalous nucleation far from equilibrium 
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We present precision Monte Carlo data and analytic arguments for an asymmetric exclusion 
process, involving two species of particles driven in opposite directions on a 2 x L lattice. We 
propose a scenario which resolves a stark discrepancy between earlier simulation data, suggesting 
the existence of an ordered phase, and an analytic conjecture according to which the system should 
revert to a disordered state in the thermodynamic limit. By analyzing the finite size effects in detail, 
we argue that the presence of a single, seemingly macroscopic, cluster is an intermediate stage of 
a complex nucleation process: In smaller systems, this cluster is destabilized while larger systems 
allow the formation of multiple clusters. Both limits lead to exponential cluster size distributions 
which are, however, controlled by very different length scales. 

PACS numbers: 64.60.Cn, 64.75.+g, 68.43.Jk 

Introduction. From recent studies of statistical model systems far from thermal equilibrium, it is clear that even 
their stationary states pose many new challenges. The well established theoretical machinery for systems in thermal 
equilibrium is essentially powerless in these contexts, and our equilibrium-trained intuition often misleads us. For 
example, though the specific heat of an equilibrium system cannot be negative, the internal energy of a non-equilibrium 
steady state may decrease when coupled to a thermal bath at a higher temperature p|. In general, non-equilibrium 
steady states depend sensitively on the details of the microscopic dynamics, resulting in widely diverse behavior at 
the macroscopic level. So far, an overarching framework is still lacking, but much insight has been gained by studying 
simple model systems. 

One such class of models are driven diffusive systems Q. Close relatives of the asymmetric exclusion process 
they involve interacting diffusing particles, driven into selected spatial directions by an external force. Models of this 
kind have been invoked to describe vehicular and pedestrian traffic gel electrophoresis [f|, and a wide range of 
biological problems, ranging from molecular motors [|| to protein synthesis 0. A particularly rich phase diagram 
emerges for two species of particles, driven in opposite directions on a two-dimensional (2D) periodic lattice Even 
with no interaction except excluded volume, we observe transitions from a freely flowing homogeneous state to a 
"jammed" state displaying a macroscopic cluster of particles. These transitions persist even if neighboring particles 
are allowed to exchange places, with a small rate 7 [s|. In contrast, an analytic solution 10] for the same system in 
ID shows that the homogeneous state prevails always, for any 7 > and particle density p < 1. Exploring notions 
of a lower critical dimension, we studied a quasi-lD system, involving just two "lanes", i.e., 2 x L, with L up to 10 4 
[Hill- Remarkably, for small 7 and half-filling, the jammed state re-emerged, with the length of the jam scaling 
with L. However, it was argued recently 0] that this jam is a mere finite-size effect and should break up into a 
disordered state when L exceeds a characteristic crossover length L c . While the latter is not known exactly, it might 
be as large as 10 70 0! 

Since simulations are an essential tool in the study of noncquilibrium steady states and should provide a fairly 
accurate picture of the thermodynamic limit, such stark discrepancies between simulation data and analytic arguments 
are disconcerting and need to be explored further. In this letter, we report extensive, high precision computer 
simulations, studying the full (p, 7) space and a wide range of L for this two-lane system. Defining an "s-cluster" as s 
particles of either species, connected to each other via nearest-neighbor bonds, we focus on the normalized residence 
distribution, p(s), which is the probability that a randomly chosen particle belongs to an s-cluster. For a homogeneous 
state, p(s) is monotonically decreasing. In contrast, for a jammed system it displays two distinct peaks: one at the 
origin (s = 1) and another marking the macroscopic cluster, at s cx L If a crossover length L c exists, this second 
peak should disappear when L increases beyond L c . Expecting L c to decrease with 7, we consider much larger 7's 
than in previous studies jll| and find, indeed, that systems with p > 0.3 cross over to monotonic distributions with 
increasing L. However, at lower densities (e.g., p — 0.1) we discover the opposite behavior: here, p(s) is monotonic 
at small L and develops a second peak as L becomes larger. In order to interpret these complex finite-size effects, we 
propose a new picture, in which larger clusters compete with a homogeneous background of "travellers" , consisting 
only of isolated particles and very small clusters. In the remainder of this letter, we briefly review the model, present 
the simulation results, and provide the analysis that leads us to our conclusions. 

The model. In a 2 x L lattice with fully periodic boundary conditions, each site can be either empty (0) or occupied 
by a single particle, which we identify as either "positive" (+) or "negative" (— ). Configurations are labeled by {a XtV }, 
where a takes on values 0, ±1 and x £ [1, £]; y = 1,2. The two species are "driven" in opposite directions, much 
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FIG. 1: The residence distribution p(s) vs. s on lattices with L = 2 10 , 2 12 , 2 14 and 2 16 for (a) p = 0.5 and 7 = 0.45 and (b) 
p — 0.1 and 7 = 0.163. The inset on the left figure magnifies the small-s data and the inset on the right shows the change from 
a double-peaked to a single-peaked distribution by increasing L. 



like cars on a 2-lane road. Starting from an initial configuration with equal numbers of the two species occupying a 
fraction, p, of the sites (i.e., N + = N- = pL), we evolve the system by exceedingly simple rules. In one Monte Carlo 
step (MCS), the following steps are repeated 2L times: (i) a pair of two nearest-neighbor sites (a "bond") is chosen at 
random; (ii) for bonds along the L direction, a (+0) or a (0—) pair is always exchanged while a (H — ) pair is exchanged 
with probability 7; (Hi) for transverse bonds, particle-hole pairs are always exchanged while particle-particle pairs are 
exchanged with probability 7. Other moves are not allowed. 

To probe the system's behavior and the effects of finite size, we explore a wide range of 7 and p, and lattices with 
L = 2 8 , 2 10 , 2 12 , 2 14 , and 2 16 ~ 64K. Since our main focus is the steady state, we discard up to 2 x 10 6 MCS 
before taking measurements. Thereafter, quantities of interest are recorded every 100 MCS, for up to 10 7 MCS. On a 
64-bit machine, we exploit a fast multispin coding algorithm (e.g., 0|) by encoding the state on a site for 64 different 
lattices (realizations) simultaneously. Thus, we are confident that our statistical errors are minimal. 

Simulation results. We set the scene by presenting a characteristic data set for the residence distribution p(s) at 
p = 0.5 and 7 = 0.45 (Fig.^,), in order to demonstrate the crossover from bimodal to single-peaked distributions with 
increasing L. The smallest system, L = 2 10 , shows a distinct peak at s ~ 380, while the next size, L — 2 12 , shows only 
a shoulder which broadens into a well-defined exponential, oc exp [— s£(7, p)}, with a 'slope' £(0.45, 0.5) ~ 0.25 x 10~ 3 
for the largest system. It is remarkable that, even for such large L's, the large s data still depend quite sensitively 
on L. On the other hand, as the inset shows, the small s part is manifestly independent of L, a feature that persists 
in the whole (7, p) space. Since we expect the crossover length L c to decrease with 7, we consider systems with 
7 = 0.5, keeping p fixed at 0.5. Here, even the smallest system (L = 2 10 ) already displays a monotonic distribution. 
In contrast, if 7 = 0.4, the distributions remain double-peaked for all L < 2 16 . In fact, for those L, the (second) peak 
position scales as s ~ 0.8 L 095 which one might (naively) interpret as the signature of a "macro scop ic" cluster, unless 
even larger L are investigated. To summarize, these findings are consistent with the conjecture |l3| ] that the presence 
of 'macroscopic' clusters is a finite-size effect. For L — > 00, we may expect residence distributions to comprise two 
components: one at small s (which we associate with the "travellers") and a slow exponential decay exp [—5^(7, p)] 
(associated with large scale "jams"). 

Surprisingly, the behavior of our model is richer than this simple picture suggests. For smaller density and 7, 
there is an additional crossover, in the opposite direction: Distributions are monotonic in small L systems and become 
bimodal as L increases, only reverting to being monotonic for much larger L. The first of these crossovers is illustrated 
in Fig. which shows residence distributions for the same four L's but at p — 0.1 and 7 = 0.163. We re-emphasize: 
In contrast to Fig. ^i, here the second peak emerges with increasing L. 

Fig- E suggests that, for any given L, it is possible to identify two distinct regions in (j,p) space: In one of them, 
located near the large 7, small p corner, p(s) is monotonically decreasing in s; in the other region (small 7, larger p), 
p (s) is bimodal. To summarize our data, the corresponding boundaries, 7 c (p, L), are drawn in Fig. [21 for several L. 
If larger systems favor (disfavor) bimodal distributions, j c (p,L) shifts up (down) with L. 

Analytic arguments. In the following, we present a simple scenario which encompasses the full range of these 
phenomena. Our discussion begins at very low densities. A space-time plot (Fig. [3]}) for p = 0.1, 7 = 0.1, and 
L = 2 8 shows the evolution of a typical configuration. With most of the system occupied by travellers, local density 
fluctuations occasionally "nucleate" a single larger cluster, containing about 24 ± 8 particles, which survives for up 
to 1000 MCS, then disintegrates and reforms elsewhere. When one of these large clusters is present, the background 
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density of travellers is depleted noticeably. Reducing 7 lengthens the cluster life times, until a single cluster persists 
for the whole run. In contrast, at higher 7, large clusters just blink on and off, and the traveller component dominates. 

Motivated by these observations, we describe the system at such low densities in terms of two "components" . One 
of these is purely uniform, consisting only of travellers at density p. The other exhibits a single large cluster, of mass 
s, reducing the density of travellers in the remainder of the system to p(s) — [2Lp — s) / (2L — s). In the spirit of 
mean-field theory, we neglect fluctuations in the cluster size as well as any charge asymmetries. Since the total number 
of each species is conserved, it is natural to focus on, e.g., the current of positive particles through either lane, as a 
function of the local densities. Of course, the current in the interior of a cluster, j c (s), should display a very different 
density dependence from the current carried by the travellers, jt (s). Assuming the traveller region to be approximately 
uniform, jt(s) should be well represented by its mean-field form [jj], jt{s) = p(s) [1 — p(s)] /2 + (7/4)p 2 (s). Here, the 
first term reflects particle-hole exchanges, proportional to the density of positive particles, p(s)/2, and that of holes, 
1 — p(s). The second term models the charge exchanges. For the current through a cluster of size s, the exact form 
is known only for the one-lane model, j c {s) = (7/4) [l + b/s + (3(s -2 )] with b — 3/2 0|. For the two-lane system, 
we have to rely on simulations which confirm the same form up to a slightly modified coefficient, b ~ 1.6 |l3l Il6|. 
Even though the 0(s ) corrections are not taken into account, this form appears to be sufficiently accurate for our 
purposes. 

Clearly, a large cluster of size s can co-exist with a uniform region only if the two currents, j c (s) and jt(s), are 
equal. If j c {s) > jt{s), the cluster will lose particles and shrink; for j c (s) < jt{s), it will grow in size. This allows 
us to identify both stable and unstable solutions of j c {s) = jt{s), as functions of 7, p, and L, in the physical domain 
s > 1. Leaving details to a future publication |17|. we just outline the main qualitative predictions here. To begin 
with, we note that j c (s) is independent of p and L while jt(s) varies only very weakly with 7. In other words, we can 
plot jt{s) for a given (p, L) and then watch j c (s) shift up or down, proportional to 7. In this fashion, we can easily 
identify three regimes: (i) At small p and L, and 7 sufficiently large, no solution exists: j t {s) remains well below j c {s), 
and the system is homogeneously filled with travellers. While local fluctuations occasionally generate small clusters, 
these disintegrate almost immediately. The two smallest system sizes of Fig. ^> belong into this regime, (ii) As 7 
decreases, the two curves jt(s) and j c (s) first touch, and then two solutions emerge: an unstable one (s u ) close to 
the origin which plays the role of a critical droplet, and a stable one at s > s u which sets the cluster size. If s is 
not much larger than s u and fluctuations are included, we expect to see purely uniform configurations compete with 
those supporting one larger cluster, as in Fig.|2K. For those parameter values, we predict s = 24, in remarkably good 
agreement with the data. Returning to Fig. we note that two solutions can also emerge if L increases at fixed 7, 
as illustrated by the L — 2 12 data. Analytically, we expect s = 70 here. Since the two currents differ by less than 5 
up to s ~ 120, clusters of all these sizes form quite easily, giving rise to the observed p(s). (Hi) Finally, as 7 is lowered 
further, only the stable solution survives and shifts to s ~ 0(L). Before this occurs, however, fluctuations begin to 
play a much more important role: Considering the L = 2 16 system in Fig. our theory would predict s a — 1436, 
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FIG. 3: Space-time plot at p = 0.1, 7 = 0.1, L = 2 8 (a) and p — 0.4, 7 = 0.4, and L — 2 9 (b). The spatial (temporal) coordinate 
runs horizontally (vertically, in units of 50 MCS). A black dot indicates a particle in either lane. 



but the largest observed clusters remain well below 700 particles. Monitoring typical configurations, we discover that 
the system is attempting to nucleate a second cluster in the traveller region, in addition to a persistent first one. 

To summarize, our simple mean-field theory describes the "nucleation" of a single cluster, from a uniform background 
of travellers, remarkably well. Naturally, once a single cluster exists, and either p or L are increased (or 7 reduced), we 
should expect, and do indeed observe, the "nucleation" of additional clusters in the traveller region. As this process 
continues, the residence distribution broadens, and eventually crosses over to a slow exponential decay, as in Fig.^. 
To illustrate this, we present a space-time plot for p — 0.4, 7 = 0.4 and L — 2 9 (Fig. [3J}) which shows several clusters 
of different sizes in each time slice. 

Conclusions. We have presented extensive simulation data for a 2 x L system with two species of particles driven in 
opposite directions. Varying L, density p, and the particle exchange rate 7, we monitor the cluster size ("residence") 
distribution p(s). We find two well-separated domains in (p, 7) space: one at low p, high 7 with monotonically 
decreasing p(s), and one at high p, low 7 with bimodal distributions. The boundary between these domains shifts 
with L, the general trend favoring monotonic p(s) as L increases. 

To understand the underlying physics at the crossovers, we draw some insight from nucleation processes. For any 
fixed value of 7, systems are uniform for sufficiently low p and L, and the associated residence distribution falls off 
rapidly. As either p or L increases, a single large cluster is nucleated. The residence distribution develops a distinct 
shoulder and, eventually, a second peak. This transition is well captured by a simple mean-field theory which balances 
the current through the cluster with that through the uniform region. As p or L increase further, a second cluster is 
nucleated, and eventually a whole population is established. As a result, the residence distribution becomes monotonic 
again, exhibiting a slow exponential decay. 

With presently accessible L's, we may only conclude that our data for 7 C (p, L) (cf. Fig. is consistent with the 
conjecture [3. At the same time, given the absence of an exact solution, we cannot rule out the possibility that 
7 C (p, L) might have a non-trivial L — > 00 limit. We do stress, however, that its non-monotonic L-dependence (at fixed 
p) appears unprecedented in equilibrium nucleation [lsf or phase separation |19| processes. Whatever the eventual 
conclusion, we believe that it is important to explore finite-size effects, since these are relevant for most applications 
of interest. For example, in the modelling of protein synthesis, system sizes ranging from IK to 10M cover the typical 
lengths of mRNA molecules. In these contexts, we feel that both crossovers may be of theoretical as well as practical 
interest. 
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